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ABSTRACT 



Context. The boundary conditions of the heliosphere are set by the ionization, density and composition of inflowing interstellar matter. 
Aims. Constraining the properties of the Local Interstellar Cloud (LIC) at the heliosphere requires radiative transfer ionization models. 
Methods. We model the background interstellar radiation field using observed stellar FUV and EUV emission and the diffuse soft 
X-ray background. We also model the emission from the boundary between the LIC and the hot Local Bubble (LB) plasma, assuming 
that the cloud is evaporating because of thermal conduction. We create a grid of models covering a plausible range of LIC and LB 
properties, and use the modeled radiation field as input to radiative transfer/thermal equilibrium calculations using the Cloudy code. 
Data from in situ observations of He , pickup ions and anomalous cosmic rays in the heliosphere, and absorption line measurements 
towards e CMa were used to constrain the input parameters. 

Results. A restricted range of assumed LIC H i column densities and LB plasma temperatures produce models that match all the 
observational constraints. The relative weakness of the constraints on N(H i) and T h contrast with the narrow limits predicted for the 
H° and electron density in the LIC at the Sun, ;i(H°) = 0.19 - 0.20 cm" 3 , and ;i c = 0.07 ± 0.01 cm 3 . Derived abundances are mostly 
typical for low density gas, with sub-solar Mg, Si and Fe, possibly subsolar O and N, and S about solar; however C is supersolar. 
Conclusions. The interstellar gas at the Sun is warm, low density, and partially ionized, with «(H) ss 0.23 - 0.27 cm -3 , T ~ 6300 K, 
A"(H + ) ~ 0.2, and X(He + ) ~ 0.4. These results appear to be robust since acceptable models are found for substantially different input 
radiation fields. Our results favor low values for the reference solar abundances for the LIC composition. 
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1. Introduction 

The solar wind evacuates a cavity in the interstellar medium 
(ISM) known as the heliosphere, from which interstellar ions 
are excluded. In contrast, neutral interstellar gas flows through 
the heliosphere until destroyed by charge-transfer with the solar 
wind and photoionization. These neutrals form the parent popu- 
lation of pickup ions (PUI) and anomalous cosmic rays (ACR) 
observed inside of the heliosphere. The properties of the sur- 
rounding interstellar medium set the boundary conditions of the 
heliosphere and determine its configuration and evolution. An 
ionization gradient is expected in the cloud feeding ISM into the 
heliosphere because of the hardness of the interstellar radiation 
field and the low opacity of the ISM dCheng & Bruhweilerll990t 
IVallered[l996t ISlavin & Frischll2002l hereafter SF02). Because 
of this ionization gradient, the densities of partially ionized 
species in the local interstellar cloud (LIC) will differ from the 
values in the circumheliosphere interstellar medium (CHISM) 
that forms the boundary conditions of the heliosphere. Hence 
we have undertaken a series of studies to determine the bound- 
ary conditions of the heliosphere based on both astronomical and 
heliospheric data. In turn our results provide tighter constraints 
on the heliosphere models used to calculate the filtration factors 



for neutrals that then permit comparisons between ISM inside 
and outside of the heliosphere. 

The distribution and velocity of interstellar H° inside 
of the heliosphere were first determined over 30 years ago 
from t he florescence of solar Lyq radiation from these 
atoms dThomas & Krassal Il97lb iBertaux & Blamontl Il97lb 



Ada ms & Frischl 1977b . Similar observations of solar 5 84 A flo- 
rescence from interstellar He showed that H/He ~ 6 for in- 
terstellar gas inside of t he heliosphere , in contrast to the cos - 
mic value H/He = 10 dAiellol 1 1978t IWeller & Meier! Il98ll) . 
More recent measurements of «(H°) compared to «(He°) in- 
side of the heliosphere find a similar ratio of H / He ~ 6 - 7 
dRichardson et all 12004 iGloeckler & Geissl 12004 Iwittel 12004 
Mobi us et al.l l2004). This difference can be attributed to two ef- 
fects: the loss of 40-60% of interstellar H° due to charge-transfer 
with protons in the helioshe ath region, a process denoted "filtra- 
tion" (Ripken & Fahd[l983T) . and the hardness of the interstellar 
radiation field at the Sun that ionizes more He than H (©. 

The Copernicus satellite first showed that the local inter- 
stellar cloud (LIC) surrounding the Sun is low de nsity, ~ 0.1 
atoms cm" 3 , partially ionized («(H + ) ~ n(H°), lYorkl 11974 
1975b, and warm (temperature < 10 4 K, e.g. 
1978). Copernicus, FUSE and HST data have 
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shown that the cluster of local interstellar clouds (CLIC, low 
density clouds within ~ 30 pc), has low column de nsities, N(H i) 
< 10 18 7 , and N ionization levels of > 30% (e.g. iLehner et al.l 
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120031: IWood et al.l20 05). indicating partially ionized gas because 
H and N ionization are coupled by charge-transfer. The low col- 
umn densities of the LIC itself, N(Hi) < 10 18 cirT 2 , indicate 
it is partially opaque to H ionizing radiation but not to He ion- 
izing radiation. Cloud opacities of unity are reached for A^(H i) 
~ 10 17 2 cirT 2 for photons close to the ionization threshold of 
hydrogen (13.6 eV), and A^(H i) ~ 10 177 cirT 2 for photons at the 
He ionization edge (24.6 eV). The result is that the LIC is par- 
tially ionized with a significant ionization gradient between the 
edge and center. Because of this radiative transfer effects are im- 
portant and need to be modeled carefully in order to determine 
the boundary conditions of the heliosphere. 

The LIC belongs to a flow of low density ISM embedded 
in the ve ry low density and apparently hot ( T ~ 10 6 K) Local 
Bubbl e dFrischl 11981b IFrisch & Yorkl 119861; iMcCammon et al l 
119831: ISnowdenet all 1 19901) . The~bulk motion of the CLIC 
through the local standard of resQ corresponds to a velocity of 
-17.0 to -19.4 km s~ ' from the direction t ~ 331°, b ~ -5° 
(Frisch & Slavin 2006). This upwind direction is near the center 
of the Loop I superbubble and the center of the "ring" shadow 
that has b een attributed to the merging of Loop I and the Local 
Bubb le (IFrisch & Yorkll986HFrischl2007l:lEpger & Aschenbachl 
1 1995b . Indiv idual cloudlets with distinct velocities are identified 
in this flow dLallement etal.lll986t IFrisch et al.ll2002l) . The ve- 
locity of the cloud feeding gas into the heliosphere has been de- 
termined by the velocity of interstellar He in the helio sphere 
measu red by the GAS detector on Ulysses, -26.3 km s" 1 dWittd 
120041) . 

Curiously, absorption lines at the LIC velocity are not ob- 
served in the nearest interstellar gas in the upwind hemisphere, 
such as towards the closest star a Cen or tow ards 36 Oph lo- 
cated ~ 5 pc beyond the heliosphere nose ([A dams & Frisch 
19771 Landsman et al.lfl984l iLinskv & Woodlll996t IWood et al. 
2000a b). This lack of an absorption component at the LIC ve- 
locity in the closest stars in the upstream direction indicates that 
the Sun is near the edge of the LIC, so the CHISM may vary 
over short distances (and hence timescales). 

We are able to test for possible past variations by compar- 
ing the first interplanetary H° Lya glow spectrum obtained by 
Copernicus in 1975 with Hubble Space Telescope observations 
of the Lya spectra obtained during the mid-1990's solar min- 
imum conditions. The observed H° velocity and intensity has 
not varied to within uncertainties over the twenty-year period 
separating these two sets of observations, so that the CHISM 
velocity field is relatively sm ooth over spatial scale s of ~ 120 
AU in the downwind direction (IFrisch & Slavinl2005l) . The 1975 
Copernicus data were acquired in the direction corresponding 
to ecliptic longitudes of A = 264.3°, = +15.0°, or ~ 13.3° 
from the most recent upwind dire ction derived from SOHO H i 
Lya data (Ouemerais et al. 2006b a, Q06). The Copernicus look 
direction was just outside of the "groove" expected in the Lya 
glow in the ecliptic during solar minimum. The groove is caused 
by increased charge-transfer in the solar wind curr ent sheet, 
which has a small tilt during minimum conditions dBzowskil 
120031) . The velocity of the Lya profile observed by Copernicus 
corresponds to -24.8 + 2.6 km s _1 , after correction to the SOHO 
upwind direction. Since Q06 measured H° upwind velocities 
during the solar minimum years of 1996 and 1997 of -25.7 + 0.2 
km s _I and -25.3 + 0.2 km s _1 , which is consistent with the 
Copernicus results, we must arrive at the conclusion that the flow 
of interstellar H° into the heliosphere was relatively constant be- 



1 Heliocentric motions are converted to the local standard of rest, 
LSR, using the Standard solar apex motion. 



tween the years of 1975 and 1997, so that any observations of 
variations in the interplanetary Lyes' glow properties must be due 
to solar activity properties alone. The thermal broadening of the 
Copernicus spectrum corresponded to a temperature of ~ 5400 
K, however measurement uncertainties allowed temperatures of 
up to 20, 000 K. 

In this paper we present new photoionization models of the 
LIC and show they reproduce both the densities of the ISM at the 
heliosphere and column densities in the LIC component towards 
the star e CMa. In earlier papers, photoionization models of both 
the "LIC" and "Blue Cloud" components obse rved towards e 
CMa were used to constrain the models (SF02, Frisch & Slavinl 
2003). SF02 grouped the properties of the LIC and Blue Cloud, 
both of which are with 2.7 pc of the Sun because they are also 
observed towards Sirius, ~ 12° from the e CMa sightline. The 
LIC velocity and density are sampled by observat ions of inter - 
stella r He i inside of the h eliosphere (-26.3 km s"',|Witte 2004), 
and iGrv & Jenkins! d2001l) show that the properties of the LIC 
and Blue Cloud differ somewhat. The present study therefore fo- 
cuses on obtaining the best model of heliosphere boundary con- 
ditions by using only data on the LIC inside of the heliosphere, 
and the LIC component towards e CMa. The present study also 
benefits from new atomic rates for the critical Mg n— >Mg i di- 
electronic recombination coefficient, improved cooling rates in 
the radiative transfer code Cloudy, recent values for the pickup 
ion densities at the termination shock, and recent values for solar 
abundances. 

The most truly unique quality of the LIC is that we are inside 
of it and therefore have the ability to sample the cloud directly 
via in situ observations carried out by a variety of spacecraft 
within the Solar System. Of all the available measurements of 
LIC gas flowing into the Solar System, the observations of the 
density and temperature of neutral He are apparently the most 
robust. Helium, unlike hydrogen, undergoes little ionization or 
heating in traversing the heliosheath regions, no deflection due 
to radiation pressure, and is destroyed by photoionization and 
electron-impact ionization within ~ 1 AU of the Sun, so we ex- 
pect that the density and temperature of He derived from the 
observations in the Solar System are tr uly representative of the 
values in the LIC (M obius et alj |2004). In this paper we put a 
special emphasis on matching these Hei data by determining 
model parameters that yield close agreement with the n(He°) and 
7\He°) data simultaneously. 

2. Photoionization Model Constraints and 
Assumptions 

The primary data constraints on our photoionization mod- 
els are the LIC com ponent column densities towards e CMa, 
dGrv & Jenkins 200jl hereafter GJ01), in situ observations 
of He (IWlttel 12004 iMobius et all I2004T). pickup ions (PUI, 
Gloeckler & Fisk 2007), and anomalous cosmic rays (ACR, 
ICummings et alj|2002l) . The astronomical and in situ observa- 
tional constraints are summarized in Table Q] 

2.1. Astronomical Constraints - The LIC towards e CMa 

The data toward e CMa of GJ01 show four separate veloc- 
ity components detected in several different ions including C n, 
Sim, Sin, and Mgn. One of these velocity components, with 
a heliocentric velocity of ~ 17 km s _1 , is identified as the LIC. 
Another is identified with a second local cloud, the "Blue Cloud" 
(BC) at ~ 10 km s . In our previous study of LIC ioniza- 
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Table 1. Observational Constraints 
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/V(Hi)//V(Hei) 


14 + 0.4 C 
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T(K) 


6300 ± 340 
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n(He°) (cm- 3 ) 


0.015 ± 0.003 




4, / = 1 


i 


n(N°) e (cm- 3 ) 


5.47 ± 1.37 x 10- 6 


3,/ 


= 0.68 - 


0.95 d 


;j(O ) e (cm" 3 ) 


4.82 ±0.53 x 10- 5 


3,/ 


= 0.64 - 


0.99 d 


;j(Ne ) e (cm- 3 ) 


5.82 ± 1.16 x 10- 6 


3,/ 


= 0.84 - 


0.95 d 


n(Ar°) c (cm- 3 ) 


1.63 + 0.73 x 10- 7 


3,/ 


= 0.53 - 


0.95 d 



a N(C ii*), N(N i), N(0 1), N(Mg n), /V(Si n), 7V(S n) and N(Fe ii) are 
used to constrain the input abundances of the models. 

b The first number indicates the reference (below), and "S" indi- 
cates that the column density is based on a saturated line. Filtration 
factors (/, 92.2| > are listed for neutrals that have been observed as 
pickup ions. 

c The uncerta i nty gi ven is only that due to uncertainties listed in 
Dupui s et al.1 dl995l) for the observed H i and He i column densities 
with the implicit assumption that the ratio is the same on all lines of 
sight. Given the substantial intrinsic variation in this ratio, however, 
the quoted uncertainty must be regarded as a lower limit to the true 
uncertainty. 

d We assume in this paper that for He / Hc = 1 , and / is not allowed 
to exceed 1. Heliosphere models predict the range /n e = 0.92 - 1 

for He and f H = 0. 50 - 0.74 for H (see filtration factors in 

Cummings etail 12002) ; iMuller & Zankll2004d : iMuller et all 120071 : 
2(3). 



Izmodenov et all 2004 
Note that the pickup ion ratios are for values at the termination 
shock. N°, O , Ne°, and Ar° densities need to be corrected for fil- 
tration in the heliosheath regions i2.2\ . 



References: (1) iGrv & Jenkins! {2001) Note t hat the values show n are 
for the LIC component tow ards e CMa.; (2) lDupuis et all ( 1 19951) : (3) 
Gloeckler & Fisld d2007l) : (4) IWittel d2004l) . 



tion (SF02) we noted that since the BC is cl ose to the LIC, 
and both are detected tow ards Sirius at 2.7 pc (Lallemen t et all 
Il994t iHebrardet ail 1 19991) . perhaps actually abutting the LIC, 
we should treat the two as a single cloud for this line of sight. 
Counter arguments to this line of reasoning include the fact that 
the BC apparently has different properties t han the LIC, though 
whether it is colder or hotter is unclear (see Heb rard et al.l l 1999, 
GJ01). For this reason and because we wish to find the ioniza- 
tion models that best predict interstellar neutral densities inside 
and outside of the heliosphere, we assume in this paper that the 
LIC and BC are truly separate clouds and select only the LIC 
components towards e CMa to constrain the models. 

In the context of these models, the best astronomical con- 
straints on the neutral ISM component are N i, the saturated O i 
line, and to some extent Mg i (though the line is weak). The elec- 
tron density can be deduced from the ratio Mg n/Mg i, which is 
determined by photoionization and both dielectronic and radia- 
tive recombination, and the excitation of the Cn fine-structure 
lines, Cii/Cii*. However, the heavy saturation of the Cn 1335 



A line in the ISM limits the accuracy of the determination of 
yV(C ii), so in this paper we have de-emphasized C n/C n* as an 
ionization diagnostic. For a discussion of the use of C n/C n* and 
oth er observations in diagnosing the C abundance and C/S ratio 
see ISlavin & Fr isch (2006). Elements with first ionization po- 
tentials (FIP) < 13.6 eV (e.g. Mg, Si, S, and Fe) are generally 
almost entirely singly ionized in the LIC and thus the column 
densities of these ions are close to the total column densities for 
the elements. 



2.2. In situ Constraints ■ 
Cosmic Rays 



He , Pickup Ions, and Anomalous 



Neutral atoms in the LIC penetrate the outer heliosphere regions, 
and become ionized primarily by charge-transfer with the so- 
lar wind ions, phot oionization, and electron impact ionization 
dRucinski et al.l 1 1 9961) . The composition of this neutral popula- 
tion reflects the partially ionized state of the LIC, rather than 
indicating a pure FIP effect. Thus the observed abundances of 
neutrals representing high FIP elements, such as He, Ne, and Ar, 
as well as H, N, and O, do not reflect their elemental abundances 
directly. As a result, these neutrals provide an interesting and 
unique constraint on the photoionization models. Once ionized, 
these interstellar wind particles form a population of ions with 
a distinct velocity distribution that are "picked up" by the solar 
wind and conve rted outwards, where they are measured by var - 
ious spacecraft dMobius et al.l l20Q4t iGloeckler & Geissl 12004). 
These pickup ions (PUI) are accelerated in the heliosheath re- 
gion and form a population of cosmic rays with an anomalous 
composition reflective of their origin as inter stellar neutrals in 
a partially ionized gas (Cummings et al. 2003). In situ observa- 
tions of these byproducts of the ISM interaction with the helio- 
sphere, H, He, N, O, Ar, and Ne, provide a unique opportunity 
to constrain theoretical models of an interstellar cloud using the 
combination of sightline-integrated data, and data from a "sin- 
gle" spatial location, the heliosphere. 

We adopt the Ulysses He measurements as the best set of 
constraints on the ISM inside of the heliosphere. The Ulysses 
satellite provides direct measurements of interstellar He at 
high e cliptic latitud es throughout the solar cycle (the GAS de- 
tector, IWittd [2004) and also measur ements of the He pickup 
ion component (the SW ICS detector IGloeckler & Geissl 12004 
IGloeckler & Geiss 2007h . Although interstellar He close to the 
Sun is detected through the resonant scattering of solar 584 A 
radiation, geocoronal contamination o f the interstellar sig nal is 
present so we prefer the Ulysses data (Mobius et al. 2004). We 
adopt the Ulysses GAS and SWICS results, n(He°) = 0.0151 + 
0.0015 cm" 3 , 7/(He°) = 6, 300 + 340 K. 

Data on the density of neutral N, O, Ne, and Ar in the sur- 
rounding ISM are provided by PUI and ACR data. The densities 
of interstellar N°, O l) , Ne°, and Ar° at the termination shock are 
listed in Table [TJ These densities must be corrected by the fil- 
tration factors, which correspond to the ratios of the densities 
at the termination shock to those in the LIC. Filtration occurs 
when neutral interstellar atoms are removed from the inflow by 
charge-transfer with interstellar protons as the atom crosses the 
heliosheath re gions. Filtration values a re listed in Table [T| base d 
on values in ICurnmings et all d2002l) : IMuller & Zaiikl d2004bl) : 
llzmodenov et alj ( 20041) . Fikration values larger than 1 are not 



considered, although some models suggest possible net creation 
of O through ch arge-transfer between Q + and H° in the he- 
liosheath regions (Mull er&Zaiikll2004al) . We adopt / He = 1 for 
He. Our models must be compared to the interstellar densities 
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log Wavelength (A) 



Fig. 1. The modeled interstellar radiation field at the Sun (model 
26) is shown as a function of wavelength (bottom X-axis) and 
energy (top X-axis). The black histogram is the modeled hot 
gas (i.e. Local Bubble) spectrum while the gray histogram is 
the cloud boundary contribution. The other line is the stellar 
EUV/FUV background. The list of elements at the top of the 
plot identifies the ionization potentials for neutrals and ions of 
interest. The energy/wavelength at which an optical depth of 1 is 
reached for several different H i column densities is shown along 
the bottom of the plot. Observed flux levels of the soft X-ray dif- 
fuse background in the Wisconsin Be and B band s are plotted as 
lines dBloch et alJll986tlMcCammon et al.lll983l) . 

obtained by correcting densities at the termination shock by fil- 
tration factors. 

2.3. Interstellar Radiation Field at the Cloud Surface 

The spectrum and flux of the cosmic radiation field control 
the ionization of the very local warm partially ionized medium 
(WPIM). The local interstellar radiation field (ISRF) that ionizes 
the LIC and other nearby cloudlets is determined by the location 
of the Sun in the interior of a hole in the neutral interstellar gas 
and dust referred to as the Local Bubble. The clustering tendency 
of hot O and B stars, and attenuation of radiation by interstellar 
dust and gas, yield the well known spatial variation of the in- 
tensity and spectrum of the ISRF. We ignore possible tem poral 
variations of the radiation field (e.g. lParravano et ai] l2003). and 
model the ISRF throughout the LIC based on the observations 
of the present-day radiation field. 

The radiation field we use in our models is based on obser- 
vations of the ISRF at the Sun, supplemented by theoretical cal- 
culations of the spectra in the EUV and soft X-rays where lack 
of sensitivity and/or spectral resolution require the use of mod- 
els to create a realistic spectrum. The directly observed radiation 
field includes the far ultraviolet (FUV) field created primarily 
by B stars, the extreme ultraviolet (EUV) field from two B stars 
(e CMa and /3 CMa) and hot nearby white-dwarf stars, and the 
diffuse soft X-ray background (SXRB). We show the radiation 
field at the cloud surface in Fig.Q] Because the interstellar opac- 
ity for radiation with E > 13.6 eV is vastly greater than for 
E < 13.6 eV, the EUV part of the ISRF originates in nearby 
regions with N(Hi) 10 18 cirT 2 while the FUV comes from a 
much larger volume. 

We use the EUV field of lVallereal dl998l) . which is based 
on data collected by the Extreme Ultraviolet Explorer (EUVE) 



satellite. The EUVE spectrometers were sensitive over the wave- 
length range of 504 - 730 A and showed that the stellar part of 
the EUV background is dominated by e CMa and /3 CMa with 
substantial contributions from nearby hot white dwarfs at shorter 
wavelengths. Vallerga extrapolated those measurements to the 
H° ionization edge at 912 A using a total interstellar H° density 
towards e CMa of JV(Hi) = 9 x 10 17 cm 2 . This value for jV(H i) 
appear s somewhat high based on observations of iGrv & Jenkins! 
(2001) which, though dependent on assumptions for gas phase 
abundances, indicate a value for A^(Hi) of ~ 7 x 10 17 cm -2 . This 
uncertainty in the total N(}\ i) affects only the extrapolated por- 
tion of the spectrum in our models since the other portion of the 
spectrum is derived by de-absorbing the observed spectrum by 
the value of N(}\ i) assumed just for the LIC. For the results pre- 
sented here we have assumed that the total N(}\ i) towards e CMa 
is 7 x 10 17 cm -2 . Assuming a larger value would increase the flux 
in the extrapolated region, though our calculations indicate that 
the overall affect on our results is only at the ~ 2% level at most. 

The FUV field is important because it sets the ionization rate 
of Mg°, which has a first ionization potential of 7.65 eV (1621 
A). The radiation field shortwards of 1600 A is heavily domi- 
nated by O and B stars in Gould's Belt, particularly those oc- 
cupying the unattenuated regions of the third and fourth galac- 
tic quadrants. A pronounced spatial asymmetry in the 1565 A 
radiation field has been observed by the TD-1 satellite S2/68 
telescope survey of the interstellar radiation field, and we use 
those data and the extrapolation down to 912 A from the 15 64 A 
measurements as calculated by Gondhalek ar et al.l (Qj80). The 
asymmetries in the TD-1 1565 A radiation field are repro duced 
by diffuse interstellar radiation field models (Henry 2002). 

The diffuse soft X-ray background (SXRB) has been ob- 
served over the entire s ky at relatively low spatial and spectral 
resolution by ROSAT ( Sno wden et al.l IT9971) and proportional 
count ers flown on sounding rockets by the Wisconsin group 
(e.g.. iMcCammon et al.l [l983l) . The broadband count rates in 
the low energy bands, particularly the B and C bands (130— 
188 eV and 160-284 eV respectively) have been modeled as 
coming from an optically thin, hot plasma at a temperature of 
~ 10 K that occupies the low density cavit y extending to ~ 50 - 
200 pc from the Sun in all directions (ISnowden et al.l [1990). 
Refinements to this picture have been required by ROSAT data 
showing absorption by relatively dist a nt clo uds (e.g. MBM12, 
ISnowden etal.lll993l) . ISnowden et all (1998) propose a picture 
in which the emission is divided between a LB component (un- 
absorbed except for the LIC) and a distant absorbed component, 
mainly in the Galactic halo. More recently there has been grow- 
ing evidence that some, possibly large, fraction of the SXRB is 
generated within the heliosphere from charge exchange between 
solar wind ions and neutral atoms (SWCX). We discuss this fur- 
ther in fl5~6l 

We model the spectrum measured by the broad-band soft X- 
ray observations by assuming that the SXR background consists 
of a local (unabsorbed) component and a distant (halo) com- 
ponent absorbed by an Hi column density of 10 19 cm -2 . The 
emission measure or intrinsic intensity of the local and distant 
componen ts are assumed to be equal . The spectrum is calculated 
using the Raymond & S mith! (1 19771 updated) plasma emission 
code assuming a hot, optically thin plasma in collisional ion- 
ization equilibrium. We explore temperatures for the hot gas of 
log T/, = 5.9, 6.0 and 6.1. The total flux, scaled by the emission 
measure, is fixed so that the B band flux matches the all-sky av- 
erage from lMcCammon et al.l d 1 983b . 
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The boundary region between the warm LIC and the adja- 
cent Local Bubble plasma may be another significant source of 
EUV radiation and we include this flux in our models 1-30. We 
model this transition region as a conductive interface between 
th e LIC and Local Bubbl e plasma in the same way as described 
in ISlavin & Frischl (120021) . The cloud is assumed to be steadily 
evaporating into the surrounding hot gas. The partially ionized 
gas of the LIC is heated and ionized as it flows into the Local 
Bubble. The ionization falls out of equilibrium, with low ion- 
ization stages persisting into the hot gas. The non-equilibrium 
ionization is followe d and the emission i n the boundary is calcu- 
lated again using the lRavmond & Smith! (I 19771) code. Ions in the 
outflow are typically excited several times before being ionized 
and the boundary region radiates strongly in the 13.6 - 54.4 eV 
band. The contribution of the interface emission to the total B 
band flux is taken into account in the calculation of the emission 
measure for the hot gas of the Local Bubble so that the B band 
flux still matches the all-sky average. 

We note that no attempt is made to make this model con- 
sistent with the size of the local cavity, which is proposed to 
contain the hot gas i n the standard model for the SXRB (e.g., 
ISnowden et all 19981) . The pressure in the hot gas is not adjusted 
to fit such models. Rather the total pressure in the hot gas for an 
evaporating cloud model is dictated by the assumed density, tem- 
perature and magnetic field in the cloud. In fact the pressures in 
our models come out far too low for the standard model to pro- 
duce the SXRB within the confine s of the Local Cavity as de- 
duced, e.g., from Na i observations (Lallement et al. 2003). This 
in turn means that if the thermal pressure in the Local Bubble 
turns out to be much lower than was assumed in those models 
because a substantial fraction of the SXRB comes from SWCX 
then the emission from the cloud boundary would be unaffected. 

An example of the profiles of the hydrodynamic variables 
for two different models is shown in Fig. [2] The magnetic field 
strength in these calculations is assumed to be proportional to 
the density at every point in the outflow. The treatment here is 
the same as in ISlavinl ([19895 which contains a more thorough 
discussion of the issues involved in this sort of calculation. The 
effect of the field on the conductivity is parametrized in a simple 
way by a constant reduction factor of 0.5. The importance of the 
magnetic field for our calculations lies in the way it affects the 
thermal pressure in the layer. Since the density drops sharply in 
the outflow and |B| ~ n, any magnetic pressure (~ B 2 ) drops 
even more sharply. The total pressure is roughly constant in the 
boundary, so the thermal pressure necessarily rises to make up 
for the decreasing magnetic pressure. Since all our models have 
nearly the same thermal pressure in the cloud, the primary effect 
of the magnetic field is to help determine the thermal pressure 
in the interface and thus radiative flux from the boundary. Since 
the total soft X-ray flux is fixed by requirin g a match with the 
Wisc onsin B band all-sky average count rate (McCammon et al. 
119831) . a larger assumed magnetic field increases only the EUV 
flux which is not constrained by the B-band data. The affect on 
the cloud of a larger EUV flux is to increase its temperature and 
ionization. 

A secondary effect of the magnetic field can be seen in the 
temperature profiles in Figure [2] In evaporating clouds in the 
ISM, if the temperature gradient is large enough and the thermal 
pressure is low enough, the conduction becomes "saturated", 
which means that the heat flux expected from the gas temper- 
ature and temperature gradient exceeds the flux that can be car- 
ried by the electrons dCowie & McKedll977l) . Saturation leads 
to a steepening of the temperature gradient and a relative slow- 
ing of the mass loss rate. In the two cases shown in Figure [2j 
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Fig. 2. Profiles of hydrodynamic variables in two different cloud 
boundary calculations. The solid lines are for model 6 which has 
n H = 0.273 cm~ 3 , log T h = 6.0, B = 2^G and N m = 4.5 x 10 17 
cm -2 , while the dashed lines are for model 8, which differs from 
model 6 only in the strength of the magnetic field, Bo = 5 fiG. 
Note that in the upper right panel, the plot of thermal pressure, 
the radial scale is much smaller in order to show the variation in 
thermal pressure. The higher magnetic pressure inside the cloud 
for model 8 leads to the higher external thermal pressure for that 
case. The temperature profile differs in the two cases because the 
degree of heat flux saturation is reduced for the higher thermal 
pressure of model 8, which in turn leads to a shallower tempera- 
ture gradient. 



the lower magnetic field case (Bq = 2 yuG) is moderately satu- 
rated (in terms of Cowie & McKee's parameter, <x = 3) while 
the other case (Bo = 5 fiG), because of the higher magnetic pres- 
sure, is less saturated. 



3. Photoionization Models 

The photoionization models of the LIC are developed follow- 
ing the same underlying procedure as in SF02, but selecting 
only the LIC absorption component towards e CMa for com- 
parison. Improvements include using updated values for in situ 
ISM observations and using a recent version of the Cloudy 
radiative transfer/th ermal equilibrium code (version 06.02.09a, 
iFerland et al.lll998l) . We run Cloudy with the assumption of a 
plane-parallel cloud geometry and format our calculated pho- 
toionizing spectrum to be used as input. The selected options 
include utilizing recent calculations for the dielectronic recom- 
bination rates for Mg + -> Mg° (e.g. lAltun et al]|2006l and the 
commands "set dielectronic recombination Badnell" and "set ra- 
diative recombination Badnell") , the assumption of constant 
pressure, and the inclusion of interstellar dust grains at 50% 
abundance compared to a standard ISM value. As we discuss 
below, the only role that dust plays is net heating of the gas since 
there is far too little column for extinction to be important (ex- 
pected E(B — V)~ 10~ 4 ' 2 ). The fraction of the heating provided 
by dust is ~ 4% of the total heating, while dust provides ~ 2% of 
the cooling (mainly by capture of electrons onto grain surfaces). 
A cosmic ray ionization rate at the default background level of 
2.5 x 10 17 s _1 (for H ionization) is included. The LIC is assumed 
to be in ionization and thermal equilibrium, and Cloudy calcu- 
lates the detailed transfer of radiation, including absorption and 
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Table 2. Model Input Parameter Values 







Input Parameter 








log T h 


So 


N m 


Model No. 


(cm- 3 ) 


(K) 


(y"G) (10 


17 cm" 2 ) 


1 


0.273 


5.9 


2.0 


3.0 


2 


0.273 


5.9 


2.0 


4.5 


3 


0.273 


5.9 


5.0 


3.0 


4 


0.273 


5.9 


5.0 


4.5 


5 


0.273 


6.0 


2.0 


3.0 


6 


0.273 


6.0 


2.0 


4.5 


7 


0.273 


6.0 


5.0 


3.0 


8 


0.273 


6.0 


5.0 


4.5 


9 


0.273 


6.1 


2.0 


3.0 


10 


0.273 


6.1 


2.0 


4.5 


11 


0.273 


6.1 


5.0 


3.0 


12 


0.273 


6.1 


5.0 


4.5 


13 


0.218 


5.9 


2.0 


3.0 


14 


0.218 


5.9 


2.0 


4.5 


15 


0.218 


5.9 


5.0 


3.0 


16 


0.218 


5.9 


5.0 


4.5 


17 


0.218 


6.0 


2.0 


3.0 


18 


0.218 


6.0 


2.0 


4.5 


19 


0.218 


6.0 


5.0 


3.0 


20 


0.218 


6.0 


5.0 


4.5 


21 


0.218 


6.1 


2.0 


3.0 


22 


0.218 


6.1 


2.0 


4.5 


23 


0.218 


6.1 


5.0 


3.0 


24 


0.218 


6.1 


5.0 


4.5 


25 


0.226 


5.9 


4.7 


3.0 


26 


0.213 


5.9 


2.5 


4.0 


27 


0.226 


6.0 


3.8 


3.0 


28 


0.216 


6.0 


2.1 


4.0 


29 


0.232 


6.1 


3.4 


3.0 


30 


0.223 


6.1 


0.05 


4.0 


42 


0.218 


6.1 




4.5 



scattering of the radiation incident on the cloud surface, as well 
as the diffuse continuum and emission lines generated within the 
cloud. 

The procedure for creating a model begins with generating 
the incident radiation field at the cloud surface ( 32.31) . The ra- 
diative transfer model is then run, and the output predictions of 
the model are compared to observations of interstellar absorp- 
tion lines in the LIC towards e CMa (7V(Cn*), N(Ni), N(Oi), 
N(Mg n), N(S n), 7V(Si n), 7V(Fe n)) and in situ observations of 
«(He°) by spacecraft inside of the solar system. The abundances 
of C, N, O, Mg, Si, S and Fe are then adjusted to be consistent 
with the observed column densities. With these new abundances, 
the Cloudy run is repeated and this process is continued until no 
adjustment of the abundances is needed. Because the abundances 
do have an impact on the emission from the cloud boundary, the 
cloud evaporation model is then re-run with the new abundances 
as well to re-generate the input ionizing spectrum. The iterative 
process of generating the spectrum and doing the Cloudy pho- 
toionization runs generally requires only a couple runs of the 
cloud evaporation program and a few runs of Cloudy. 



4. Model Results 

To investigate the dependence of the results on the input pa- 
rameters we calculate a grid of 24 models. We explore total H 
density, «h = 0.273,0.218; Local Bubble hot gas temperature, 
log Ti, = 5.9, 6.0, 6.1; magnetic field strength, Bo = 2, 5 fiG; and 
Hi column density, N m = 3 x 10 I7 ,4.5 x 10 17 cirr 2 . The val- 
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Fig. 3. Model results for the He density and temperature in the 
ISM just outside the heliosphere. The squares, circles, and trian- 
gles are for models that are part of the initial grid of 24 models, 
while the stars are for models 25 - 30 for which the magnetic 
field, B, and the total H density, n H , were varied to match the ob- 
served n(He°) and T. For the grid models, the empty symbols are 
for models with N(Hi)= 3 x 10 17 crrr 2 , while the filled symbol 
models have N(Hi)= 4.5 x 10 17 crrr 2 . As the legend shows, the 
color (black vs. gray) indicates the magnetic field strength and 
the symbol shape indicates the temperature assumed for the hot 
gas of the Local Bubble. For identical kinds of points, the one to 
the left is for a model with «(H°) = 0.218 crrr 3 and the one to 
the right has n(H°) = 0.273 cm 3 . The ellipse is the error range 
around the observed values for n(He°) and T. 



ues for T and n(He i) at the solar location, the endpoint of the 
calculations, are shown in Fig.[3]for this set of models. We then 
explore another six models in which we do the calculations over 
a grid in log 7), and Nhi but vary the values of «h and Bo in or- 
der to match the observed T and n(He i). We employed a multiple 
linear regression to assist in narrowing down the search region 
for the values of «h and Bo needed to match the observations. 
Since the dependencies of the results on the parameters really 
are quite non-linear, this procedure could not work to predict 
exactly correct values for the required parameters, but was use- 
ful for getting close to the correct values. Based on the results 
for the initial grid of models, we use log 7), = 5.9,6.0,6.1 and 
N m = 3 x 10 17 ,4 x 10 17 crrr 2 for this smaller grid, models 
25 - 30. We chose to use 4 x 10 17 crrr 2 rather than 4.5 x 10 17 
crrr 2 because the higher column density models, for the most 
part, produced temperatures that were too high. As can be seen 
from Fig. [3] all of these models (25-30), plotted as stars, are con- 
sistent with the observed T and n(He i), indicated by the ellipse 
in the figure. Table|2]gives the input parameters for each model. 
Model predictions for the e CMa sightline integrated through the 
LIC are presented in Table [3] Model predictions for the CHISM 
(i.e. at the solar location) are shown in Table[4] 

In SF02 we tested models with no emission from the cloud 
boundary. This amounts to assuming that the boundary is a sharp 
transition from the hot gas of the Local Bubble to the warm gas 
of the LIC. We have again explored such models using the LIC 
data as constraints in the same way as for the models discussed 
above with a conductive interface. When there is no evaporative 
boundary, our models do not depend on the magnetic field in the 
cloud since in that case the ionizing flux consists only of dif- 
fuse emission from the hot gas of the Local Bubble and EUV 
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Model" 


log/yCHtot) 


log N(At i) 


log A'(Arn) 


logA'(Sim) 


A<(Mg 11) 
iV(Ma 1) 


N(C 11) 

N(C II* ) 


A I H 1 ) 
W(Hc I) 


Obs. b 


- 


- 


- 


12.40 


^"-i in 


93 - 430 c 


12 - 16 


1 


17.58 


11.48 


11.76 


8.895 


706.3 


213.0 


11.36 


2 


17.78 


11.67 


11.95 


9.641 


325.3 


234.1 


11.40 


3 


17.65 


11.33 


11.84 


9.844 


246.1 


182.1 


11.44 


4 


17.82 


11.55 


12.00 


10.12 


136.5 


198.0 


11.80 


5 


17.58 


11.49 


11.75 


8.961 


717.8 


230.0 


11.62 


6 


17.78 


11.67 


11.94 


9.711 


286.7 


244.9 


11.62 


7 


17.64 


11.34 


11.81 


9.909 


207.4 


194.8 


12.47 


8 


17.82 


11.56 


11.97 


10.17 


119.7 


208.2 


12.66 


9 


17.59 


11.46 


11.75 


9.289 


553.0 


239.5 


12.27 


10 


17.78 


11.66 


11.93 


9.861 


210.3 


248.1 


12.14 


11 


17.64 


11.33 


11.78 


9.993 


171.2 


201.7 


13.44 


12 


17.82 


11.54 


11.95 


10.23 


104.3 


210.6 


13.52 


13 


17.60 


11.44 


11.78 


9.100 


768.9 


255.1 


11.52 


14 


17.79 


11.64 


11.96 


9.738 


336.4 


274.4 


11.47 


15 


17.67 


11.30 


11.85 


9.944 


251.2 


215.9 


11.65 


16 


17.84 


11.51 


12.01 


10.21 


140.6 


232.1 


12.01 


17 


17.58 


11.47 


11.76 


8.926 


840.1 


262.5 


11.65 


18 


17.79 


11.65 


11.95 


9.792 


308.9 


286.3 


11.64 


19 


17.66 


11.30 


11.82 


10.02 


208.5 


229.2 


12.76 


20 


17.84 


11.52 


11.98 


10.26 


123.5 


242.1 


12.95 


21 


17.59 


11.43 


11.76 


9.318 


633.1 


280.8 


12.41 


22 


17.79 


11.63 


11.94 


9.926 


224.9 


288.4 


12.25 


23 


17.66 


11.29 


11.79 


10.10 


175.0 


235.5 


13.81 


24 


17.84 


11.50 


11.96 


10.31 


110.2 


245.5 


13.85 


25 


17.66 


11.32 


11.84 


9.887 


271.3 


213.6 


11.71 


26 


17.74 


11.56 


11.91 


9.676 


383.9 


270.1 


11.70 


27 


17.63 


11.36 


11.79 


9.751 


331.1 


242.0 


12.51 


28 


17.73 


11.59 


11.90 


9.625 


411.6 


284.5 


11.79 


29 


17.62 


11.38 


11.77 


9.720 


329.3 


251.6 


13.01 


30 


17.72 


11.60 


11.88 


9.642 


386.6 


295.1 


12.08 


42 


17.77 


11.70 


11.92 


9.555 


482.9 


317.9 


11.51 



a The best models, consiste nt with all observationa l data (see jj4), are indicated by bold face. 

b Observational results from Gry & Jenkins (2001) (see tablefj}. The value s listed for W(Hi)/A'(Hei) are the range of values observed excluding 
Feige 24 which is one of the most distant stars observed bv lDupuis et alj i 19951) and has unusually large i V(H i) and ratio values. 

c The upper limit on A'(Cii) is not well determined observationally because of the saturation of the line. iGrv & Jenkins! (|2001) define it by 
assuming a solar C/S abundance ratio and using the observed S n column density. We find (see lSlavin & Frischl2006l) that the abundance of C 
required to match the N(C if ) observations is supersolar, with C/S~ 36 - 44, which results in a much higher upper limit on N(C n). 



emission from nearby hot stars and the spectrum is not related 
to the properties of the cloud. For these models our model grid 
consists of total H density, «h = 0.273, 0.218, Local Bubble hot 
gas temperature, logJ/, = 5.9,6.0,6.1, and Hi column density, 
Nm = 3xl0 17 ,4.5xl0 17 cm~ 2 , for a total of twelve models that 
we label as models 31^-2. Of these models, however, only two 
resulted in ionizing fluxes sufficient to heat the cloud to a tem- 
perature ~ 6000 K at the same time as matching the constraints 
on the ion column densities. These models were the ones with 
logT h = 6.1 and N m = 4.5 x 10 17 cnr 2 (models 36 and 42). 
For the other cases, either at the surface or deeper into the cloud, 
there are insufficient photons to provide the heating to balance 
the cooling and the cloud temperature drops sharply to < 1000 
K. The model with «h = 0.218 (model 42) is consistent with the 
observed value of n(He°) and with the column density data as 
well. 

From these 42 models we have selected those that provide 
acceptable results for the observational constraints, according to 
prioritized requirements. The first requirement is that the model 
predict the density and temperature of He i observed inside of 
the solar system. Models 14 (marginally), 15 and 25-30 and 42 
predict a He density and temperature consistent with the ob- 



served values within the reported errors. They also predict the 
PUI Ne densities, for an assumed Ne/H=123 ppm. Models 26 
and 28 successfully match the PUI data for Ne/H as low as 
~ 100 ppm. These models are also required to match the ob- 
served Mgn/Mgi ratio in the LIC towards e CMa. Models 14, 
27, and 29 marginally fit this criterion, while 26, 28, 30, and 42 
successfully fit this criterion. Note that the new models now pro- 
vide acceptable predictions for the CHISM temperature, which 
was not the case for the best models in SF02. These models are 
also consistent with the observed C n/C n* ratios in the LIC com- 
ponent towards e CMa though this constraint is weak because of 
the large uncertainties in N(C n). Based on these comparisons, 
and given the uncertainties in both data and models, Models 14, 
26-30, and 42 are plausible models, but Models 26 and 28 appear 
to best match the observational constraints. 

Based on the constraints and assumptions presented in the 
previous section, we select models 14, 26-30 and 42 as the best 
models for the LIC ionization, with models 26 and 28 favored 
by the PUI Ne data provided that Ne/H > 100 ppm. We believe 
that these seven models then give a realistic range for the uncer- 
tainties in the boundary conditions of the heliosphere, providing 
that the underlying assumptions implicit in the Cloudy code, e.g. 
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Table 4. Model Results for Solar Location 
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a The best models (see © a re indicated by bold face. 

b Observational results from Gloeckler & Geiss (2004), Gloeckler (2005, private communication) and Witte et al. ( 1996) (see table[TJfor uncer- 
tainties). 



photoionization equilibrium, are correct. The predictions of the 
best models provide excellent matches to the observational con- 
straints. Based on these models, we find the boundary conditions 
of the heliosphere to be describable as n(H°) = 0.19 -0.20 cm -3 , 
« e = 0.05-0.08 cm" 3 , andX(H) = H + /(H° + H + ) = 0.19-0.26, 
X{He) = He + /(He° + He + ) = 0.36 - 0.43. For these models, we 
find abundances of O/H = 295 - 437 ppm, C/H = 589 - 813 
ppm, and N/H = 40.7 - 64.6 ppm (Table©. The total LIC den- 
sity is n(H)o = 0.213 - 0.232 cm- 3 , while the strength of the 
interstellar magnetic in the cloud varies between and 3.8 fiG. 
The Ne PUI data further favor densities of n(H°) « 0.19 cm" 3 
and « e = 0.06 - 0.07 cm 3 . 

In this analysis we have assumed negligible filtration for He 
in the heliosheath regions. Modeling of the He filtration factor 
however allows values as small as /n e = 0.92, which yields 
w(Hei)= 0.0164 cm 3 for the CHISM (Table Q). The predictions 
of Model 21 agree with this value for n(Hei), as well as with 
the ratios Mg 11/Mg 1 and C 11/C 11* towards e CMa and the pickup 
ion Ne and Ar data (Table |3)- The predicted cloud temperature 
is low by ~ 1000 K. The density and ionization for Model 21 are 
n(H°) = 0.21 cm 3 and n e = 0.06 cm 3 . We therefore conclude 
that our best models listed above are robust in the sense that they 
predict consistent values for the n(H°) to within 5%, and electron 
densities to within 25%. 



The radiation field incident on the LIC for Model 26 is shown 
in Fig. [TJ and the spectral characteristics of the field for each 
model are listed in Table [5] The wavelengths regions A < 912 A 
and A ~ 1500 A are of primary importance for the photoioniza- 
tion of the cloud, the former because it determines H° and He 
ionization, and the latter because it determines Mg° ionization. 

The ionization parameter is defined as U = 0/n(H)c, where 
<1> is the H ionizing photon flux, n(H) is the total (neutral + ion- 
ized) hydrogen density at the cloud surface and c is the speed 
of light. The total ionizing photon fluxes at the cloud surface 
(photons cm -2 s" 1 ) for the three bands 13.6-24.6, 24.6-54.4, and 
54.4-100 eV, are given by Oh, 0^0, and <t>He+, respectively. The 
ratio of the total number of H° and He ionizing photons in the 
incident radiation field is given by 2(He°)/2(H°). The quantity 
(E) (eV) is the mean energy of an ionizing photon, equal to the 
integrated energy flux from 13.6 to 100 eV divided by the inte- 
grated photon flux over the same energy range. 

In Figure[3]we show n(He°) and temperature of the CHISM 
for our model calculations. In Figure [4] we show Nii/Ni vs. 
Mgn/Mgi, illustrating an anti-correlation of the ratios caused 
by the fact that Mg 11/Mg 1 decreases with electron density while 
N 11/N 1 indicates the ionization level in the cloud. The ionization 
of the CHISM is listed in Table for Model 26, where com- 
monly observed elements are listed. The abundances of He, Ne, 
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Table 5. Characteristics of the Model Radiation Field 
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Fig. 4. Model results for A^(N n)/iV(N i) versus Mg n/Mg i. The 
symbols have the same meaning as in Figure [3] In this case the 
models with higher n(H°) lie below and slightly to the left of 
those with lower density. A^(N ii)/N(N i) is an indicator of cloud 
ionization fraction while Mg n/Mg i goes as l/n e . The dotted line 
is the observed value for Mg n/Mg i and the dashed lines indi- 
cate the 1 — cr error range for the value. We see that models that 
match the observed ratio all correspond to relatively low ioniza- 
tion, X(H) ~ 0.20 - 0.27 in the CHISM. 



Na, Al, P, At, and Ca were assumed, based on solar abundances, 
and were not adjusted in the modeling process. The abundances 
of C, N, O, Mg, Si, S and Fe were adjusted for each model to 
match observed column densities towards e CMa (see ^2). The 
elemental abundances that have been assumed, with the excep- 
tion of that for He, are not expected to have any significant im- 
pact on the model results. 

5. Discussion 

5.1. Heliosphere Boundary Conditions 

As discussed above, the best models of those we calculated 
are determined by the match to the CHISM He density and 
temperature found by the in situ Ulysses measurements ( 92.21 . 
combined with the matching the LIC component column den- 
sity ratios Mgn/Mgi and Cii/Cii*. These models span a fairly 
large range in the model parameters: «h = 0.213 - 0.232 
cm" 3 , logT/, = 5.9 - 6.1 K, B = 0.05 - 3.8 yuG, and ^(Hi) 
= 3.0 x 10 17 - 4.5 x 10 17 cirr 2 . Despite this, the predicted val- 
ues for neutral H density and electron density in the CHISM lie 
within a remarkably small range: n(H°) = 0.19 - 0.20 cm -3 , 
« e = 0.05 - 0.08 cirr 3 , for models that include the conductive 
boundary. For the one model without evaporation that is consis- 
tent with the data we find n(H°) = 0.186 cm 3 , n e = 0.052 cm 3 . 
Including the PUI Ne data as a constraint narrows the density re- 
sults to n(H°) a 0.19 cm 3 and n e = 0.06 - 0.07 cm 3 . For these 



10 Slavin and Frisch: Boundary Conditions of the Heliosphere 

Table 6. Model 26 Results for Ionization Fractions 3 Table 7. Elemental Gas Phase Abundances (ppm) 



Element 


PPM 


I 


II 


III 


IV 


H 


10" 


0.776 


0.224 


— 


— 


He 


10 5 


0.611 


0.385 


4.36(-3) 


— 


C 


661 


2.68(-4) 


0.975 


0.0244 


0.000 


N 


46.8 


0.720 


0.280 


8.52(-5) 


0.000 


O 


331 


0.814 


0.186 


4.71 (-5) 


0.000 


Ne 


123 


0.196 


0.652 


0.152 


2.79(-6) 


Na 


2.04 


1.47(-3) 


0.843 


0.155 


6.34(-6) 


Mg 


o.ol 


1 OQ^ 1\ 

l.yo(-i) 


U.53U 


0.146 


U.UUU 


Al 


0.0794 


5.37(-5) 


0.976 


0.0118 


0.0123 


Si 


8.13 


4.21(-5) 


0.999 


8.02(-4) 


3.10(-5) 


P 


0.219 


1.35(-4) 


0.977 


0.0232 


9.29(-5) 


S 


15.8 


6.47(-5) 


0.971 


0.0288 


1.95(-6) 


Ar 


2.82 


0.263 


0.500 


0.238 


2.83(-6) 


Ca 


4.07(-4) 


9.21(-6) 


0.0155 


0.984 


1.87(-4) 


Fe 


2.51 


7.01(-5) 


0.975 


0.0245 


5.75(-6) 



11 Numbers less than 10 3 are written as x(y) where y is the exponent 
and x is the mantissa (or significand). 



densities to stray out of this range would appear to require sig- 
nificant errors in the underlying comparison data, e.g. n(He°) in 
the CHISM, or substantial non-equilibrium ionization effects in 
the LIC. The variation in the interstellar radiation field between 
the different models that match the data give us some degree of 
confidence that these densities are not highly sensitive to the de- 
tails of the radiation field. The range « e = 0.05 - 0.08 cm -3 
corresponds to an electron plasma frequency of 2.0-2.5 kHz, 
which is the frequency of the mysterious weak radio emission 
detected beyond the termination shock in th e outer heliosphere 
dGurnett & Kurthll2lMlMitchell et aljllool . 

5.2. Hydrogen Filtration Factor 

Tracers of H° inside of the termination shock, after filtration, in- 
clude the H i Lya backscattered radiation, H pickup ions, and 
the slowdown of the solar wind at distances beyond 5 AU from 
mass-loading by H PUIs. The range of n(H°) found above to 
best fit the combined heliospheric n(He°) and LIC data towards 
e CMa, «(H°) = 0.19 - 0.20 cm 3 , represents the density of neu- 
tral interstellar H atoms outside of the heliosphere, and removed 
from heliospheric influences. The hydrogen filtration factor, fu, 
can be obtained from comparisons between these models and in- 
terstellar H° densities at the termination shock as inferred from 
in situ observations of interstellar H inside of the heliosphere. 

The accompanying papers in this special section provide es- 
timates of the interstellar H° density at the termination shock. 
The solar wind slows down due to massloading by interstel- 
lar H, yielding w(H°)= 0.09 + 0.01 cirT 3 at the termination 
shock ( Richardson et al. 2007). The density of H pickup ions ob- 
served by Ulysses is inferred at the termination shock, yielding 
n(H°)= 0.1 1 + 0.01 cirT 3 ; models of H atoms traversing the he- 
liosheath regions then yield for the CHISM n(H° )= 0.20 + 0.02 
cm" 3 and» D = 0.04+0.02 cm 4 , orn e ~ 0.05 crrT 3 (Bzowski et al. 
120071) . The radial variation in the response of the interplanetary 
Lya 1215 A backscattered radiation to the solar rotational mod- 
ulation of the Lya "beam" that excites the florescence yields 
»(H°)~ 0.085-0.0 95 cm -3 , depending on the heliosphere model 
( Prv or et all 120071) . From these n(H°) values at the termination 
shock, we estimate that 43%-58% of the H-atoms successfully 
travers e the heliosheath region, or /h ~ 0.43 - 0.58. lMuller et alj 
(2007) evaluate filtration using five different plasma-neutral 



Element 



Model No. 
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Mg 


Si 
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Fe 


14 


589 


40.7 


295 


5.89 


7.24 


14.1 


2.24 


21 


955 


60.3 


447 


9.77 


11.5 


22.9 


3.55 


25 


631 


66.1 


437 


7.76 


10.0 


19.5 


3.09 


26 


661 


46.8 


331 


6.61 


8.13 


15.8 


2.51 


27 


759 


64.6 


437 


8.71 


10.7 


20.9 


3.31 


28 


708 
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331 


7.08 


8.32 


16.6 


2.57 


29 


813 


64.6 
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9.33 


11.0 


21.9 


3.39 


30 


741 


46.8 


331 


7.41 


8.51 


17.0 


2.63 


42 


724 


39.8 


295 


6.76 


7.76 


15.1 


2.34 



models, and find a range of /h = 0.52 - 0.74. A hydrogen filtra- 
tion of fu = 0.55 + 0.03 is consistent with both in situ data and 
radiative transfer models. 

5.3. Gas-Phase Abundances 

The LIC photoionization models are forced to match the ob- 
served set of column densities (Tables Q] and 0. The gas-phase 
abundances of most elements are treated as free parameters that 
can be varied in order to match observed column densities, so 
that the successful models yield elemental abundances for the 
LIC that are automatically corrected for unobserved H + (Table 
17). The exceptions are that He, Ne, and Ar abundances, being 
unconstrained by the observations toward e CMa, are not ad- 
justed but are assumed to be 10 5 ppm, 123 ppm and 2.82 ppm, 
respectively. In the models, N(Cif) is a constraint on both the 
C abundance (in place of the heavily saturated Cn 1335 A line) 
and « e , such that the product of the abundance and « e is more 
tightly limited than either quantity individually. The requirement 
to match both A^(C n*) and n(He°) effectively restricts the ioniza- 
tion fraction of H, which in turn limits O and N ionization whose 
ionization fractions are tied by charge-transfer to the H ioniza- 
tion at LIC temperatures. 

Early studies showed that the abundances of refractory el- 
ements in the very local ISM are enhanced compared to abun- 
dances in col d disk gas (Marschall & Hobbs 1972; Stokes 1978; 
iFrischl 119811) . Throughout warm and cold disk gas, the un- 
derabundances of refractory elements compared to solar abun- 
dances (by factors of 10~' - 10~ 4 ) a re taken to represent deple - 
tion onto interstellar dust grains (e.g. Savage & Sembach 1996). 
This view is supported by the correlation found between elemen- 
tal depletions and the temperature ch aracteristic of condensa- 
tion at solar pressure and composition (Ebel 2000), and assumes 
that there is a reference abundance pattern that characterizes the 
cloud, and remains constant over the cloud lifetime as atoms are 
exchanged between the gas and dust phases. Below Q5.41 i we 
compare solar abundances with observed gas-phase abundances 
to predict the gas-to-dust mass ratios for the LIC, based on the 
assumption that LIC gas and dust have remained coupled over 
the cloud lifetime. 

An important question is whether the LIC has solar abun- 
dances. Isotopes of is O and 22 Ne isotopes measured in the ACR 
population suggest that this is so. ACRs are characterized by 
a rising particle flux for energies below 10-50 MeV/nucleon, 
and this characteristic spectral signature is seen for lfi O, ls O, 
20 Ne, and 22 Ne. Ratios of 16 0/ 18 ~ 500 and 20 Ne/ 22 Ne ~ 
13.7 are found for both the ACRs and solar material, indicat- 
ing that the CHISM and solar material have similar composi- 
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Table 8. Solar Abundances (ppm) 





Grevesse 


Holweger 


Lodders 


Grevesse 




Sauval 






et al. 




(1998) 


(2001) 


(2003) 


(2007)" 


c 


334 ± 46 


391(+110,-86) 


290 ± 27 


275 ± 34 


N 


84 ± 12 


oJ.j |9 


82 ± 20 


67 ± 10 


O 


683 ± 94 


545 +l07 

-SQ 


579 ± 66 


513 + 63 


Nc 


121 ± 17 


ioo + |^ 


91+21 


77 ± 12 


Mg 
Si 


38 ±4 
35 + 4 


34.5 + ' J 
34.4+H 


41 ±2 
41 ±2 


38 + 9 
36 + 4 


S 


22 ±5 




18 ±2 


15 + 2 


Ar 


2.54 ± 0.35 




4.24 ± 0.77 


1.70 + 0.34 


Fe 


32 ±4 


28.1- 


35 ±2 


32 + 4 



a Protosolar abundances are obtained by increasing the photospheric 
ab undances by 0.05 dex for elements heavier than He, as suggested 
bv lGrevesse et al] J2007h 



tions dLeske et al.ll2000l: lLeskell200(jh . We therefore adopt solar 
abundances as the underlying reference abundance pattern for 
the LIC. 

Unfortunately a prominent uncertainty exists in the correct 
solar abundances of volatile elements such as O and S, which 
have low condensation temperatures (r con d, 180 K and 700 K 
respectively), and noble elements such as Ne and Ar. Solar 
abundances are determined from photospheric data (C, N, O, 
Mg, Si, S, Fe), the solar wind (Ar, Ne), solar active regions 
(Ne), and helioseismic data (He); abundance s of non-volatile el- 
emen t s are also found f r om meteoritic data (|Grev esse & Sail vail 
1998; iHolweeerl 120011: iLoddersi 12003b IGrevesse et aiT l2007). 
Solar abundances from these studies are listed in Table [8] Our 
results, discussed below, indicate that if the LIC has a solar abun- 
dance composition, as indicated by the 18 Q and 22 Ne data, then 
the lower abundances found by IGrevesse et al] d2007l) are pre- 
ferred by our models. 

Ne: In th ese models we have assu med the Ne abundance 
is 123 ppm dAnders & Grevesse] 1 1989b . which is based on a 
combination of photospheric and interstellar data. Solar sys- 
tem Ne abundances are difficult to meas ure because of FIP 
effects, however values include 77 ppm dGrevesse et al.| [2007. 
after adding 0.05 dex to account for gravitational settling of 
the elements), and ~ 41 ppm for solar wind in coronal holes 
dGloeckler & Geissll20"07l) . The predicted densities of n(Ne°) in 
the CHISM for models 26 and 28 (Tabled are in agreement with 
the most recent PUI results for Ne (TableQ]), and Ne densities as 
low as ~ 100 ppm are allowed when filtration is included. 

The CHISM Ne abundances indicated by these results ap- 
pear to be consistent with Ne abundances and ionization levels 
in the global ISM. The Ne a bundance in the Orion nebula is 100 
ppm dSimpson et alj|2004l) . iTakei et all d2002l) measured X-ray 
absorption edges formed by Ne and O in the interstellar gas and 
dust towards Cyg X-2, and found abundances of Ne/H ~ 92 ppm 
and O/H ~ 579 ppm when bo th atomic and co mpound forms 
in the sightline were included. lJuett et all £2006) observed the 
X-ray absorption edges of Ne and O towards nine X-ray bina- 
ries which sampled both neutral and ionized warm material, and 
found that the ionized states formed in the ionized material have 
the ratio Neiii/Nen ~ 0.23. This value is identical to the pre- 
dictions of model 26, Nem/Nen ~ 0.23, a fortuitous agreement 
that may indicate that the EUV radiation field in the CHISM is 
similar to the generic galactic EUV field in the solar vicinity. 



Ar: Solar Ar abundance determinations range between Ar/H 
~ 1.4 - 5.0 ppm (Table[8j; we have assumed Ar/H = 2.82 ppm. 
The predicted Ar density at the Sun is within the uncertainties 
of the PUI data, although the range of possible filtration factors 
(0.64 - 0.95, see 92. 2\ also allow considerable leeway. 

O: In the warm ISM such as the LIC, the ionization of oxy- 
gen and hydrogen is tigh tly coupled over timesca les of ~100 
years by charge transfer dField & Steigmanlll97ll) . so that the 
assumed A^(H i) combined with A^(0 1) measurements act to con- 
strain the deduced O abundance in the gas. The two best models 
(26, 28) correspond to O/H = 331 ppm, however the O column 
density measurements are based on the saturated 1302 A line, and 
have ~ 35% uncertainty. The modeled LIC value of O/H ~ 331 
ppm indicates that ~ 35% of the O atoms are depleted onto dust 
grains. An oxygen filtration factor of fo ~ 0.75 is required by 
the PUI data and Model 26. 

These models yield gas-phase O abundances that are con- 
sistent with observations of more distant interstellar sightlines. 
The ratio jV( Q i)/jV(H i) is measure d in both low and high extinc- 
tion clouds. lOliveira et al .1 (|2003) used unsaturated Oi lines in 
the 910-1 100A interval and found Oi/Hi = 317 ± 19 ppm for 
~ 30 sightlines that included both types of material. Sightlines 
with detected H 2 , N(U) > IP 205 cm' 2 and (» H ) = 0.1 -3.3 cm- 3 , 
yield O/H =319 + 14 ppm dMever et al.lll997t). A survey of 19 
stars with an average distance of 2.6 kpc by I Andre et alT f2003) 
found O i/H i = 408 + 14 ppm, where the long sightline and high 
average value N(H)/E(B — V) = 6.3 x 10 21 cirT 2 mag -1 indicate 
a bias towards sightlines containing many clouds that individu- 
ally have low extinctions. A larger sample of 56 sightlines for 
a range of extinctions and distances show that sightlines with 
higher average mean densities, («h), show O/H = 284 + 12 ppm, 
versus O/H = 390 ±1 ppm for stars with low values of (rin) 
(ICartledge et al.ll2004l) . For comparison, solar abundance stud- 
ies yield a range of ~ 450 - 780 ppm. 

C: The best Models 26 and 28 yield a gas-phase abundance 
of C of C/H = 661 and 708 compared to solar abundances 
of ~ 240 - 500 ppm, whic h is consistent with our earlier re- 
sults (Slavin & Frisch 2006) indicating an overabundance of C 
in the LIC. We speculate that shock destruction of carbonaceous 
grains, perhaps combined with some local spatial decoupling be- 
tween carbonaceous and silicate grains, may explain these find- 
ings. 

Singly-ionized carbon is an important coolant in the LIC 
( 95. 5b , so the C overabundance is required to maintain the tem- 
perature of the CHISM at the observed value. The carbon abun- 
dance obtained here indirectly depends on the Mg n— »Mg i di- 
electronic recombination coefficient that determines the ratio 
Mg n/Mg i, since that ratio is used as a criteria for the best mod- 
els. The same ionization correction that gives the C abundance 
also successfully predicts Ne ionization in global WPIM and the 
S abundances in the LIC, although this may be a fortuitous co- 
incidence. In the adjacent si ghtline towards Siri us the LIC has 
N(C n)//Y(H i) = 1 , 050 ppm faebrard et al.lll999b . An ionization 
correction of 300% is required to make this value consistent with 
solar abundances, and such a large ionization correction is not 
consistent with the ionization levels of X(H) ~ 20 - 26% found 
here. In contrast, sightlines wi th cold ISM show C abundances 
on th e order of 135 + 46 ppm (ISofia et al.ll 19971; ISembach et al .1 
120001) . 

N: The best models (26 and 28) find N/H = 46^17 ppm, com- 
pared to solar values of ~ 57-1 10 ppm. These results are consis- 
tent with the PUI results, N/H ~ 19-47 ppm, after filtration fac- 
tor uncertainties are included. The N and O results and favor an 
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ISM a bundance pattern for volatiles similar to the lGrevesse et al.l 
(2007) photospheric abundances. 

S: The best models predict S/H = 16-17, compared to solar 
values of ~ 13-27 (including uncertainties, see Table[8]l. Sulfur 
is found to have little or no deplet ion onto dust grains in warm 
diffuse ISM (e.g. lWeltv et"aflll999l) . 

Mg, Si, Fe: These refractory elements are observed in the 
LIC gas with abundances far below solar (factors of 3-15). 
Approximately 92%, 82%, and 77% of the Fe, Mg, and Si, re- 
spect ively, are presumab ly depleted onto interstellar dust grains. 
If dGrevesse et al J 12007) abundances are assumed for the LIC, 
then to within the uncertainties the LIC dust has the relative com- 
position of Fe:Mg:Si:0 = 1:1:1:4, as is consistent with amor- 
phous olivines MgFeSiCU. Fe and Si are dominantly singly ion- 
ized, while Mg has a significant fraction (~ 15%) that is twice 
ionized. The gas-phase abundances of these refractory elements 
are highly subsolar, even after ionization corrections are made, 
indicating that these elements are substantially depleted onto in- 
terstellar dust grains. In contrast to C, however, the silicate dust 
in the LIC that carries the missing Mg, Si, and Fe has experi- 
enced far less destruction than the carbonaceous grains. 

Can, Nai: Weak lines of the trace ionization species Can 
and Nai are common diagnostics of ionization and abundance 
for interstellar clouds, including the partially ionized LIC; Nai 
is also frequently used as a diagnostic of the H column density. 
We note that our models show that the ratios A^(Ca n)/Af(Na i), 
JV(Nai)/iV(H), and 7Y(Na i)/iV(H i) vary by 30%, 77%, and 93%, 
respectively, between the best models (Models 26-30, 42). As 
trace ionization species, the densities of Na i and Ca n are highly 
sensitive to volume density, n(Na°), n(Ca + ) oc n(H)n e . We there- 
fore conclude that Nai and Can are imprecise diagnostics of 
ionization levels, H density, and abundances in warm partially 
ionized clouds. 

5.4. Gas-to-Dust Mass Ratio 

Because the abundances are automatically corrected for unob- 
served H + , we use the model results to infer the total mass of the 
interstellar dust, providing that the gas and dust in the LIC form 
a coupled and closed system that evolves together as the cloud 
moves through the LSR. The LIC LSR velocity is 16-21 pc/Myr, 
so that a LIC origin related to the Loop I or Scorpius-Centaurus 
superbubble would require that the LIC gas an d dust remained 
a clo s ed system o ver timescales of 4-5 Myr dFrisch & Slavin! 
20061 lFrischlll981l) . Gas-to-dust mass ratios calculated from the 
best models (26 and 28) using the missing-mass argument are 
in the range Rg/d= 149 - 217, depending on solar abundances. 
The detailed information about Rg/d for different assumptions 
and the different models is listed in table [9] 

For comparison, Rg/d determined from comparisons of in 
situ observations of interstellar dust inside of the solar system, 
compared to th e gas densities of these models, yield Rg/d = 1 15— 
125 (Table l9l lLandgraf et alJl2000tlAltobelli etaIH200l . The in 
situ Rg/d is an upper limit, since the smallest interstellar dust 
grains (radii < 0.15 fim) with large charge-to-mass ratios (and 
thus small Larmor radii) are excluded from the heliosphere by 
the interstellar magnetic field which is draped over the helio- 
sphere. 

For all the models, the Rg/d determined from comparing 
in situ dust measurements with the CHISM gas mass flux is 



2 This argument assumes that the ISM reference abundances, in this 
case solar abun dances, rep resent the sum of the atoms in the gas plus 
the dust jFrisch et al. II 19991) . 



Table 9. Gas-to-Dust Mass Ratios from Models and In Situ 
Observations 



ModeP 



Source 


14 


26 


27 


28 


29 


30 


42 


GS98 


137 


149 


196 


149 


197 


150 


138 


Lodders 


158 


174 


238 


174 


239 


175 


160 


Grevesse 


194 


217 


321 


217 


323 


218 


196 


In Situ 


115 


116 


123 


116 


125 


116 


107 



The source of the comparison solar abunda nces is listed in column 
1 (Table [8] The in situ dust flux is from lLandgraf et alj (1 2000). 
correc ted downwards by 20% as recommended by Altobe lli et al.1 
d2004l) to account for side-wall impacts. 



Table 10. Major Heat Sources in LIC Gas a 



Source b 


Fraction of Heating 


Hi 


0.657 


He i 


0.248 


dust 


0.055 


Hen 


0.016 


cosmic rays 


0.010 



a Results for model 26. Other models are qualitatively the same, 
though there are some quantitative variations. 

b For lines with ion names, the source here denotes the ion that is 
photoionized. Dust heating comes from photoelectric ejection by 
photons of the background FUV radiation field. Cosmic ray heating 
comes from electron impact ionization of the gas and direct heating 
of the electrons in the LIC plasma by the cosmic ray electrons. 



lower than that determined by assuming solar abundances and 
using the gas phase abundances we determine to find Rg/d- 
This suggests that somehow the dust flowing into the helio- 
sphere is concentrated relative to the gas, compared to the over- 
all T^rc_S2gMine towards e CMa. The lower solar abundances 
of lGrevesse et al.l ([2007) result in lower required depletions, and 
produces stronger disagreements with Rg/d determined from in 
situ data. W e do not understan d this result, which we have found 
previously dFrisch et al.lll999 |). Since Rg/d is sensit ive to the 
mass of Fe in the dust grains dFrisch & SIavinl f2003). we sug- 
gest that this difference may indicate inhomogeneous mixing of 
the gas and silicate dust over the ~ 0.64 pc extent of the LIC. 

5.5. Heating and Cooling Rates 

The heating and cooling rates for Model 26 are listed in Tables 
[TOl and [TTJ The primary heat sources are photoelectrons from 
the ionization of H° and He , with dust and cosmic ray heating 
contributing less than 7% of the heating. The dominant source 
of cooling is the [Cn] 157.6 fim fine-structure line, making up 
43% of the total. This is more than twice the contribution of 
any other coolant. Nearly all the cooling is due to optical and 
infrared forbidden lines with many lines contributing at the ~ 
1 % level. H recombination, free-free emission and dust, through 
the capture of electrons onto grain surfaces, also contribute at 
about a 2% level. The importance of Cn as both a constraint 
on the C abundance in our models as well as a major coolant 
means that any model that aims to reduce the abundance of C 
to a solar level faces severe difficulties. The models with LIC 
temperatures in the T Ue o = 6 300 + 340 K range indicated by 
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Table 11. Major Coolants in LIC Gas a 



Ion/Line 


Fraction of Cooling 


[Cn] 157.6yum 


0.428 


[Sn] 6731 A 


0.145 


Fe n (total) 


0.074 


[Sin] 34.8 fan 


0.065 


[Nen] 12.8^m 


0.035 


[Oi] 63.2/jm 


0.028 


I-l rpcntnii 
STL ICCUII1U. 


n CO A 


dust 


0.024 


[Nem] 15.6ycan 


0.020 


[Nn] 6584 A 


0.018 


net free-free 


0.018 


[Oi] 6300 A 


0.017 


[On] 3727 A 


0.011 


[Am] 6.98 yum 


0.011 



a Results for model 26. Other models show similiar results. 



thermal and magnetic pressures are equal, this indicates a mag- 
netic field strength B ~ 2.7 yuG, in agreement with field strengths 
for these models. As noted is 92.31 the main effect of the field 
strength in the models is to regulate the pressure in the evapo- 
rative cloud boundary, which in turn affects the flux of diffuse 
EUV radiation incident on the cloud. The amount of EUV flux 
helps determine the temperature in the cloud, which is how the 
observational constraints fix the magnetic field strength in the 
context of our modeling. Thus we do not explicitly fix the mag- 
netic field strength with the goal of achieving equipartition and 
indeed some of our successful models have lower or higher field 
strengths. It is probably coincidental that the field strength re- 
quired to match the in situ He temperature for our best models 
is also close to the equipartition field strength, but it is at least 
encouraging that this field strength is consistent with our pho- 
toionization models. We note that if thermal, cosmic ray, and 
magnetic pressures are approximately equal the LIC has a pres- 
sure of ~ 6300 crrT 3 K. 



the in situ He data all require supersolar abundances of C. The 
total heating/cooling rate for the LIC at the Sun for this model is 
3.55 x 10~ 26 ergs cm" 3 s" 1 . 



5.6. Radiation Field 

Recently it has been proposed that a significant portion 
of the SXRB can be attributed to charge-transfer (a.k.a. 
charge exchange) between the solar wind ion s (e.g., Q +7 an d 



Q +8 ) and interstellar neu t rals dCravensI 12 000; Snowden et al 



2004; Wargelin et al. 2004; Smith et al. 2005; Koutroumpa et al 



2006). While it seems at present that some fraction of the low 
energy X-rays are from this mechanism, it is unclear how large 
that fraction is. We note that basing the properties of the local 
hot plasma in the galactic pla ne on SXRB emission at ene rgies 
E > 0.3 keV is problematical. dBellm & VaiHancourtl2 005) have 
compared the Wisconsin B and Be band data with the ROSAT 
R12 data, and concluded that the observed anticorrelation be- 
tween R12 and N(R i) indicates that more than 34% of the SXRB 
generated in the Galactic disk must come from the Local Bubble. 
They also concluded that a heavily depleted plasma with lo g 
T ~ 5.8 is consisten t both with the iMcCammon et alj (|2002); 
ISanders et al.1 (l200lh X -ray spectra l data, and the uppe r limits 
set on the EUV emission by CHIP S dHurwitz et al.l2004l) . When 
the [Robertson & Cravensl d2003l) models of SXRB production 
by charge-transfer with the solar wind are considered, then only 
half of the SXRB in the plane is required to arise from a hot local 
plasma. We also note that the atomic physics for the calculation 
of the low energy part of the emission is still quite uncertain 
(V. Kharchenko, private communication). At this point we take 
the simple approach of ignoring the charge-transfer emission, 
though we plan to consider its possible impact in future work by 
reducing the assumed SXRB flux from hot gas. As noted previ- 
ously, a lower SXRB flux due to a lower pressure in the hot gas 
does not necessarily have any impact on our calculated flux from 
the evaporative cloud boundary. 

5.7. LIC Pressure 

The strength of the interstellar magnetic field in the LIC is un- 
known, though modeling of its effects on the heliosphere sets 
some constraints. Our best models (26 and 28) presented here 
have a thermal pressure of ~ 2100 cm" 3 K for the LIC. If the 



1986 



1992; 



1999; 



5.8. Comparisons with other LISM Sightlines 

There have been a number of efforts to under stand the 
LISM ionization an d ab und ances dFrisch et al.l 
Cheng & Bruhweilerl 1 1 990t iLallement & Bertinl 
Vallergalll99d iLallement & Ferletl Il997t iHolberg et alj 
Kimura et al.l |2003j). The studies that attempt to derive gas 
phase elemental abundances find a range of results, generally 
fairly consistent with ours. A point of particular interest is the 
abundance of carbon t hat is surprisingly o verabundant in our 
results. As an example. iKimura et al.l ((2003) find (based on four 
sightlines and excluding the e CMa and B CMa sightlines), a 
subsolar C abundance in contr adiction with our re s ults. R esults 
for thirteen sightlines from iRedfield & Linskvl d2004l) with 
velocity components consistent with the LIC velocity vector 
show that N(Cv)/N(Oi) > 1 for 8 of them, especially those 
at lower column density indicating a solar or supersolar C 
abundance. For our best models N(Cii)/N(Oi) ~ 2. Our series 
of studies are unique in that we model the radiation field incident 
on the cloud, include radiative transfer effects, and calculate 
the thermal equilibrium within the cloud. The ionization varies 
through the cloud as does the temperature and density (slightly) 
and we compare observations within the heliosphere with the 
physical conditions at that point in the cloud rather than basing 
the model on line-of-sight averages. 

Our present results indicate that n(Nn)/n(Ni)~ 0.32 - 0.50 
at the solar location, with N becoming more ionized as the 
sightline approaches the cloud surface. The column density ra- 
tio is thus higher, ranging from 0.38 - 0.62. Observed values 
of N(N n)/jV(N i) toward o ther nearby stars are 0.58+"^ toward 



Cape lla (IWood et al.l2002l) . 1 .29+0.23 t oward HZ43dRruk et al.l 
2002), L91+°g towards WD1634-573 dLehner et al.ll2003l) . and 
1.13 + 0.24 towards 77 UMa (Frisch et al. in preparation). The 
total Hi column density towards each of these stars is greater 
than the N{Hi)~ 4 x 10 17 cm -2 found for the best models here. 
The nearest of these stars, Capella, has an ionization compara- 
ble to that of the LIC. The two high-latitude stars HZ43 and r\ 
UMa appear to sample low opacity regions where the ionization 
is larger than at the Sun, as does the WD 1634-573 sightline that 
appe ars to cros s the nearby diffuse Hn region seen towards A 
Sco dYorklll983l) . As we have noted the e CMa line of sight is 
special because that star is the dominant source of stellar EUV 
photons for the LIC. Thus for sightlines at a large angle from 
the e CMa sightline, if the H 1 column between points along the 
sightline and e CMa is small the apparently high column points 
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are subject to a strong EUV field. Such geometry dependent ion- 
ization effects can be important for non-spherical clouds subject 
to a strongly spatially variable ionizing radiation source. 

The variation in the fractional ionization of the CLIC gas 
has a direct impact on our understanding of the distribution and 
physical properties of low column density clouds for several rea- 
sons. (1) Abundances of elements with FIP < 13.6 eV must al- 
ways be calculated with respect to A^(Hi)+ A^(Hn) for very low 
column density clouds. (2) Cloud geometry affects the opacity 
of observed sightlines so that the opacity to ionizing radiation 
is not directly traced by the observed value of N(Hi). For lines 
of sight other than that towards e CMa, this could require more 
complex radiative transfer models in which the difference be- 
tween the line of sight toward the star and that toward one of the 
primary sources of ionizing flux, e CMa, at each point is taken 
into account. 



6. Conclusions 

There are many uncertainties regarding the detailed properties of 
the ionizing interstellar radiation field incident on the LIC. The 
data we have on the LIC, both from absorption line studies and 
in situ measurements by spacecraft in the heliosphere, provide us 
with strong constraints on the ionization and composition of the 
LIC and particularly the CHISM. By exploring a range of models 
for the ISRF we find that while a fairly broad range of radiation 
fields can produce photoionization consistent with the data, other 
outputs from the models fall within a relatively narrow range of 
values. Our results for the models explored in this paper in which 
we require our models to be consistent with the LIC component 
of the absorption lines observed towards e CMa include: 

1 . For a range of assumptions regarding the H i column density 
of the LIC, N(H i) = (3.0-4.5) x 10 17 cirr 2 , and temperature 
of the hot gas of the Local Bubble, log 7), = 5.9,6.0 and 
6.1, we are able to find model parameters that allow a match 
of the model results with best observed quantities, n(He°), 
r(He°) and N(Mgii)/N(Mgi). For these models we assume 
that the cloud is evaporating because of thermal conduction 
between the hot Local Bubble gas and the warm LIC gas and 
include the emission from the cloud boundary. 

2. For the best models in terms of fits to data, the required input 
parameters are: initial (i.e. at the outer edge of the cloud) 
total H density, n(H) « 0.21 -0.23 cirT 3 ; and cloud magnetic 
field, B < 3.8 /iG. 

3. If we assume that the magnetic field configuration reduces 
thermal conductivity at the boundary enough to prevent 
evaporation and ignore any radiation from the cloud bound- 
ary, we find that for most cases the radiation field does not 
cause sufficient heating to maintain the LIC at the temper- 
ature observed, T — 6, 300 + 340 K. One set of parameter 
choices, though, yields a successful model. These parame- 
ters are JV(Hi) = 4.5 x 10 17 cirT 2 and log T h = 6.1. 

4. Despite the wide range of possible input parameters, the out- 
put values for quantities important for shaping the helio- 
sphere are confined to a fairly small range: n(H°) = 0.19 — 
0.20 cirT 3 , and n e = 0.05 - 0.08 cm 4 . 

5. AH filtration factor of fu = 0.55 + 0.03 yields good agree- 
ment between the radiative transfer model predictions for 
«(H i) in the CHISM, and n(H i) at the termination shock as 
found from observations of PUIs, the H i Lya glow, and the 
solar wind slow-down in the outer heliosphere. This filtra- 
tion value is also consistent with heliosphere models of the 



ionization of interstellar H atoms traversing the heliosheath 
regions. 

6. Elements with ionization potentials 13.6 - 25 eV, e.g. H, He, 
N, O, Ne, and Ar, are partially ionized with ionization frac- 
tions of ~ 0.2 - 0.7. 

7. By requiring that the models match the column densities de- 
rived from absorption line data we are able to determine the 
necessary elemental abundances for several elements. We 
find that the abundances of N and O may be somewhat sub- 
solar. Sulfur is roughly solar, and C is substantially super- 
solar. Mg, Si and Fe are all sub-solar by factors of 3 - 15. 
The depletions of Fe, Mg, Si and O in the LIC are consis- 
tent with a dust population consisting of amorphous silicate 
olivines MgFeSi04, though other compositions for the dust 
are possible as well. We conclude that any carbonaceous dust 
in the LIC must have been destroyed, while silicate dust has 
persisted. Except for the gas-to-dust mass ratio, these results 
are in better agreem ent with the lower solar abundances of 
iGrevesse et alj (2007). H owev er we note that the O and Ne 
abundances o f lLoddersI (120031) are in better agreement with 
other astronomical data such as the X-ray absorption edges. 

8. The gas-to-dust mass ratio derived from missing mass in the 
gas-phase for our best models depends strongly on the as- 
sumed reference abundance set and range from 137 - 323. 
Our two best models, nos. 26 and 28, give a range of 
149-217. For these same models Rg/d = 1 15-125 based on 
the observed flux of dust into the heliosphere. The discrep- 
ancy of these values is minimized, in fact leading to con- 
sistency within the errors, if one assumes an abundance set 
such as that of GS98 which has large abundances of the met- 
als. The GS98 abundances lead to substantial O depletion, 
however, which is not easily explained and conflict with the 
S abundances found for models 26 and 28. 

9. These models also show that the densities of the trace ion- 
ization species Ca n and Na i are extremely sensitive to den- 
sity and ionization. Therefore the ratios A^(Ca n)/iV(Na i), 
N(Na i)/N(H i), and N(Na i)/N(H) are, by themselves, inade- 
quate diagnostics of warm low density diffuse gas. 
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